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abstract 

It is shown that in supersymmetric models (SUSYMs) the large supersymmetric contributions to 
B — ► X s jj, + i2^ come from the Feynman diagrams which consist of exchanging neutral Higgs bosons 
(NHBs) and the chargino-stop loop and are proportional to m^m^tai^ (3 / m\ when tan/3 is large 
and the mass of the lightest neutral Higgs boson m^ is not too large (say, less than 150 Gev). 
Numerical results show that the branching ratios of B — > X s fi + fi~ can be enhanced by more than 
100% compared to the standard model (SM) and the backward-forward asymmetry of lepton is 
o ■ significantly different from that in SM when tan/3 > 30. 

PACS numbers: 11.30.Pb, 13.20. Jf 

It is widely believed that supersymmetry (SUSY) is one of the most promising candidates for 
physics beyond SM since it offers a scheme to embed the SM in a more fundamental theory in 
which many theoretical problems such as gauge hierarchy, origin of mass and Yukawa couplings 
can be answered. One direct way to search for SUSY is to discover SUSY particles at colliders. 
But, unfortunately, so far no SUSY particles have been found. Another way is to search for 
its effect through indirect methods. In most of SUSYMs R parity is conserved so that SUSY 
contributions to an observable appear at the loop level. Therefore, it has been realized for a 
long time that rare processes can be used as a good probe for searches of SUSY, since in these 
processes the contributions of SUSY and SM arise at the same order in perturbation theory. The 
B — > X s l + l~ (l=e, [A, r) process, one of rare processes, in SUSYMs has been extensively studied 
[1-5]. The effects of large tan/3 have been noticed in recent papers [f|, ||. There is a stop-chargino 
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loop diagram which gives a large contributions to CV when tan/3 is large ||. This leads to that in 
the minimal supergravity model (mSUGRA) there are regions in the parameter space where the 
branching ratio of b — ■> sl + l~ (l=e, fi) is enhanced by about 50% compared to the SM ||. However, 
the contributions from exchanging NHBs are ignored in these previous analyses. Recently, the 
contributions of NHBs in SUSYMs have been taken into account. Because the contributions to 
b — ► st + t~ coming from the chargino-stop loop diagram are proportional to m^m T tan 3 /?/m^ 
(h = h°, A ) when tan/3 is large, the branching ratio of b — > st + t~ can be enhanced by about 
200% compared to the SM 0. 

From experimental points of view, the observation of B — > X s l + l~ (l=e, jj) is easier accessible 
than that of B — > X s t + t~ . The inclusive decay B — > has been observed by CLEO. Meantime, 
experiments at e + e~ and hadron colliders are closing in on the observation of B — > K*l + l~ (l=e, 
A*) [§]■ The B factories presently under construction will collect some 10 7 -10 8 B mesons per year 
which can be used to obtain good precision on low branching fraction modes. Therefore, it is 
meaningful to pay attention to the process B — > X s l + l~ for (l=e, fi). As pointed above, the 
contributions of NHBs are proportional to the mass of a lepton and tan 3 /3. For B — > X s e + e~, the 
contributions can be safely neglected due to the smallness of m e , no matter how large tan/3 is (of 
course, in the theoretically allowed range, say, in SUSYGUT, tan/? < 50). However, for b — > sfi + fi~, 

tan/3 can be as large as m T as long as tan/? > 17. Thus one can expect that for b — > sfi + fi~, in 
addition to the enhancement coming from the possible change of the sign of CV ( the value of Cj is 
fixed by the measurement of b — > S7 with a branching ratio of (2.32 ±0.57 ±0.35) x 10~ 4 and 95% 
C.L. bounds of 1 x 10~ 4 < Br(B — > X s 7) < 4.2 x 10~ 4 0), a even more significant enhancement 
coming from exchanging NHBs arises in SUSYMs with large tan/3. In the letter we calculate the 
invariant mass distribution and backward-forward asymmetry of dilepton angular distribution for 
B — > X s fi + fi~ in SUSYMs. Our results show that the branching ratio of B — > X s fi + fi~ is enhanced 
by at least 100% compared to the SM and the back-forward asymmetry is more sensitive to tan/3 
than the invariant mass distribution when tan/3 is larger than 30 and masses of Higgs bosons, stops 
and charginos are in the reasonable range (i.e., all constraints from phenomenology are satisfied). 
Note that the invariant mass distribution of B — > X s {i + fi~ in a two Higgs doublet model with large 
tan/? is not enhenced compared to the SM. Therefore, the rare process B — > X s fi + fi~ provides 
a good opportunity to distinguish SUSYMs with large tan/3 from the SM and it is possible that 
the first distinct signals of SUSY could come from deviations from the SM in the inclusive decay 
B -> Xsii+iir. 

Inclusive decay rates of B — > X s fi + fi~ can be calculated in the 1/mq expansion and it has 
been shown that the leading order term turns to be the decay of a free b quark and corrections 
stem from the 1 / m 2 Q order and are small about a few percent [|TIJ . Therefore in what follows we 
limit our analyses to the leading order term, i.e. the decay b — > sji + n~ . 

There are several classes of new contributions in SUSYMs and the dominant ones are arising 
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from chargino(x)-uptype squark loop and charged Higgs boson (if ± )-uptype quark loop. Because 
of small generation mixing of squarks coming from phenomenological constraints on K° — K° and 
D° — D°, the contributions from gluino-downtype squark loop and neutralino-downtype squark 
loop are much smaller than the dominant ones and are neglected in the following. 
The effective Hamiltonian relevant to the b — > sfi + fi~ process is 

A(~l 10 10 

V 1 i=l i=i 



where Oi (i=l, 2, 10) are given in Ref. [pj| , and QiS come from exchanging neutral Higgs bosons 
and have been given in Ref. [|1^]. The coefficients Ci(m w ) in SUSYMs have been calculated [[I], ||]. 
We calculate the coefficients CQi(mw) in SUSYMs and the results are: 

C Q i{m w ) = ™ bm 2n tg 2 /j{(sin 2 a + /icos 2 a)[-^(/i(x Hf ) - 
4m ft0 sm vw %wt 
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ii'=l 

where the definitions of the functions fi, Ti (i=l, 2), A±, B±, A and the meaning of the matrices U, 
V, T have been given in and we have omitted less important terms because they are numerically 
negligible compared to those given in eq. (2) when tan/5 > 20. 

From eq. (2), we see that, for large tan/3(> 20), the coefficients Cq { are proportional to 
mbm^/m,hta.n 3 (3 (h = h°,A°). One factor of tan/3 comes from the chargino-up-type squark loop 
and tan 2 /? from exchanging the neutral Higgs bosons (note that in the large tan/3 approximation 
cos _1 /3 ~tan/3). Therefore, Cq 1 can compete with Cj and even overwhelm Ci as long as tan/3 
is large enough. We remark that the chirality structure of the Qi (i=l, 2) operators allows a 
large tan/3 enhancement for the Cq 4 (i=l, 2) coefficients, as happened for the magnetic moment 
operator 7 , and there is no such a large tan/? enhancement for the Cj (i=8, 9) coefficients due to 
the different chirality structure of the (i=8, 9) operators. Incorporating the QCD corrections 
to the coefficients Cj and Cq i in the standard way, we calculate these coefficients at //=m&. 
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The differential branching ratio and the forward-backward asymmetry of the dimuon angular 
distribution for B — > X s fi + fi~ can be obtained from ]7) with substituting for m T . The numerical 
results of the invariant mass distribution and backward-forward asymmetry are shown in Fig. 1 
(a) for the mSUGRA model and Fig. 1 (b) for the MSSM with a typical choice of masses of 
sparticles and Higgs bosons respectively. The mSUGRA parameters (m , mi/ 2) A) = (190, 190, 
380) Gev, Higgs mass mixing parameter \i < and tan/?=30 have been chosen in Fig. 1 (a). In 
the computations of sparticle mass spectra and mixings we neglect the Yukawa couplings of the 
first two generations. The chosen values of masses of relevant sparticles and Higgs bosons in 
MSSM are given in the Figure Captions. The constraint from the LEP and b —>■ sj has been 
imposed in our numerical calculations. One can see from the Fig. 1 that a large enhancement of 
the differential branching ratio dT/ds shows up and the enhancement can reach 100% compared 
to SM when tan/3=30. The backward-forward is significantly different from that in SM. The 
predictions without including the contributions of exchanging NHBs are also shown in Fig. 1 in 
order to compare. It is evident from the figure that the contributions of exchanging NHBs to 
the differential branching ratio are the same order of magnitude as supersymmetric contributions 
without including exchanging NHBs in the low s region (s<0.4) and larger than those in the high s 
region (s>0.4). There are regions in the parameter space where the contributions of NHBs alone 
make a large enhancement of the differential branching ratio. For example, for a set of values of 
parameters (6^-=— 20°, m X2 =220 Gev, m Xl =100 Gev, m,j=430 Gev, 772^=250 Gev, mj 2 =500 Gev, 
77^40=80 Gev, m^ u =160 Gev, and tan/3=30) the enhancement of dT/ds coming from NHBs is 
about 80% compared to SM. 

We would like to make some remarks: 

(i) The large enhancement of the invariant mass distribution of B — > X s fi + fi~ compared to the 
SM is of a common feature of SUSYMs with large tan/? in some region of the parameter space. The 
Fig. 2 shows the results for tan/3=30. For larger tan/3, for example, tan/3=45, the enhancement 
can reach 200%. The enhancement exists as long as the mass splitting of stops is large enough 
(say, >100 Gev). The condition is necessary because if all the squark masses are degenerate 
(m^ = m~ t2 = m), the large contributions arising from the chargino-squark loop exactly cancel 
due to the GIM mechanism ||. In order to illustrate this point we show the Cq. as a function of 
the stop mixing angle 9i under the above condition in Fig. 2. As can be seen from the Figure, Cq 4 
is large enough to make an large enhancement of dT/ds in a wide range of Of (about from to 

.E) 

(ii) The QCD corrections to coefficients Cj and CQ i are incorporated in the leading logarithmic 
approximation in our numerical computations. The one-loop mixing of Qi with Or has been 
analyzed |T2| and leads to about 5% correction to Cj(mb) when tan/3=30. A next-to-leading 
order(NLO) analysis without including Qi for B — > X s l + l~ has been performed, where it is 
stressed that a scheme independent result can only be obtained by including the LO and NLO 
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corrections to Cf while retaining only the LO corrections in the remaining Wilson coefficient 



13| . Because we did not include the NLO corrections the theoretical uncertainty due to the 
renormalization scale /i dependence is about ±20% as fi is varied in the range 1/2 mb</j,<2mi ) . 
We expected that a full NLO analysis including Qi for B-^X s l + l~ will appear in the near future. 
As pointed in the Ref. citebsg5, there is a SUSY high scale uncertainty and it is possible that the 
scale /i dependence is large enough to effectively encompass the uncertainty 

(iii) The following values of parameters have been used in the numerical calculations: m t = 
175Gev, m c /mb = 0.3, r]=a s (mb) / a s {m w ) = 0.548. We have estimated the uncertainties from the 
parameters and results are that the mt dependence is weak and the uncertainties are about ten 
percent. The error from neglecting the strange quark mass m s is of order m 2 s jm\ and consequently 
is very small. 

In summary, we have investigated the differential branching ratio and backward-forward asym- 
metry of lepton for B^X s fi + fi~ in SUSYMs with large tan/3. There is a 100% enhancement of the 
differential branching ratio compared to SM if tan/?>30 and the masses of Higgs bosons, squarks 
and charginos are in the reasonable range. Because there is almost no enhancement till tan/?=50 
in a two Higgs doublet, one can make the conclusion that the first distinct signals of SUSY could 
come from the observation of B — > X s fi + {i~. 
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Figure Captions 

Fig. 1 dT/ds and A(s) for the case /i < and tan/3=30, a) = m = 190 Gev, A = 380 

GeV in the mSUGRA and b) %=-40°, m X2 =200 Gev, m Xl =90 Gev, m Q -=350 Gev, m fl =220 Gev, 
m t - 2 =450 Gev, m A o=80GeV, 77^=160 Gev in the MSSM. The solid , dashed and dotted lines 
represent the predictions of the SUSYMs, the SUSYMs without including contributions of NHBs 
and SM respectively. 

Fig. 2 Cqi and Cq2 varying with the mass splitting of stop, the mass splitting of chargino and 
the stop mixing angle 9i for /j<0, tan/9=30, 772^=150 Gev, m Xl =90 Gev and m^o = 80 Gev; the 
characters following the line style indicate the mass splittings: the first character means the mass 
splitting of stop (h represents 300 Gev and 1 100 Gev), and the second character means that of 
chargino (h represents 410 Gev, m 210 Gev, and 1 110 Gev). 
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(a)tan/3=30, m =m 1/2 =190 GeV, A=380 Gev 
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Fig. 2 
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